Non-linear QCD dynamics in two-photon interactions at high energies 



V.P. Gongalves 1 , M.S. Kugeratski 2 , E.R. Cazaroto 3 , F. Carvalho 4 and F.S. Navarra 3 
1 Institute) de Fisica e Matemdtica, Umversidade Federal de Pelotas 
Caixa Postal 354, 96010-900, Pelotas, RS, Brazil. 
2 Centro de Engenharia da Mobilidade, Umversidade Federal de Santa Catarina, 
Campus Universitdrio, Bairro Bom Retiro 89219-905, Joinville, SC, Brazil. 
3 Instituto de Fisica, Universidade de Sao Paulo, CP. 66318, 05315-970 Sao Paulo, SP, Brazil 
4 Departamento de Ciencias Exatas e da Terra, Universidade Federal de Sao Paulo, 
Campus Diadema, Rua Prof. Artur Riedel, 275, 
Jd. Eldorado, 09972-270, Diadema, SP, Brazil 

Perturbative QCD predicts that the growth of the gluon density at high energies should saturate, 
forming a Color Glass Condensate (CGC), which is described in mean field approximation by the 
Balitsky-Kovchegov (BK) equation. In this paper we study the 77 interactions at high energies 
and estimate the main observables which will be probed at future linear colliders using the color 
dipole picture. We discuss in detail the dipole - dipole cross section and propose a new relation 
between this quantity and the dipole scattering amplitude. The total 77, 7*7* cross-sections and 
the real photon structure function F^XjQ 2 ) are calculated using the recent solution of the BK 
equation with running coupling constant and the predictions are compared with those obtained 
using phenomenological models for the dipole-dipole cross section and scattering amplitude. We 
demonstrate that these models are able to describe the LEP data at high energies, but predict a 
very different behavior for the observables at higher energies. Therefore we conclude that the study 
of 77 interactions can be useful to constrain the QCD dynamics. 
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I. INTRODUCTION 



The high energy limit of perturbative QCD is characterized by a center-of-mass energy which is much larger than 
the hard scales present in the problem. The simplest process where this limit can be studied is the high energy 
scattering between two heavy quark-antiquark states, i.e. the onium-onium scattering. For a sufficiently heavy onium 
state, high energy scattering is a perturbative process since the onium radius gives the essential scale at which the 
running coupling a s is evaluated. In the dipole picture [l|, the heavy quark-antiquark pair and the soft gluons in 
the limit of large number of colors N c are viewed as a collection of color dipoles. In this case, the cross section can 
be understood as a product of the number of dipoles in one onium state, the number of dipoles in the other onium 
state and the basic cross section for dipole-dipole scattering. At leading order (LO), the cross section grows rapidly 
with the energy (a oc a 2 s e^ aip -^ Y , where (a P - 1) = ^JY^ m2 ~ q.5 and Y = In s/Q 2 ) because the LO BFKL 
equation Q predicts that the number of dipoles in the light cone wave function grows rapidly with the energy. Several 
shortcomings are present in this calculation. Firstly, in the leading order calculation the energy scale is arbitrary, 
which implies that the absolute value of the total cross section is therefore not predictable. Secondly, a s is not running 
at LO BFKL. Finally, the power growth with energy violates s-channel unitarity at large energies. Consequently, new 
physical effects should modify the LO BFKL equation at very large s, making the resulting amplitude unitary. 

A theoretical possibility to modify this behavior in a way consistent with the unitarity is the idea of parton 
saturation, where non-linear effects associated to high parton density are taken into account. The basic idea is that 
when the parton density increases (and the scattering amplitude tends to the unitarity limit), the linear description 
present in the BFKL equation breaks down and one enters the saturation regime. In this regime, the growth of the 
parton distribution is expected to saturate, forming a Color Glass Condensate (CGC), whose evolution with energy 
is described by an infinite hierarchy of coupled equations for the correlators of Wilson lines (For recent reviews see 
Q). In the mean field approximation, the first equation of this hierarchy decouples and boils down to a single non- 
linear integro-differential equation: the Balitsky-Kovchegov (BK) equation. In the last years the next-to-leading order 
corrections to the BK equation were calculated through the ressumation of a s Nf contributions to all orders, 

where Nf is the number of flavors. Such calculation allows one to estimate the soft gluon emission and running 
coupling corrections to the evolution kernel. The authors have found that the dominant contributions come from 
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the running coupling corrections, which allow us to determine the scale of the running coupling in the kernel. The 
solution of the improved BK equation was studied in detail in Ref. Q . In Q a global analysis of the small x data for 
the proton structure function using the improved BK equation was performed (See also Ref. @). In contrast to the 
BK equation at leading logarithmic a s ln(l/a;) approximation, which fails to describe the HERA data, the inclusion 
of running coupling effects in the evolution renders the BK equation compatible with them (See also |9h111]]). 

A reaction which is analogous to the process of scattering of two onia discussed above is the off-shell photon scattering 
at high energy in e + e~ colliders, where the photons are produced from the lepton beams by bremsstrahlung (For 
a review see, e.g., Ref. fl^). In these two-photon reactions, the photon virtualities can be made large enough to 
ensure the applicability of perturbative methods or can be varied in order to test the transition between the soft and 
hard regimes of the QCD dynamics. From the point of view of the BFKL approach, there are several calculations 
using the leading logarithmic approximation [l3l [l4| and considering some of the next-to-leading corrections to the 
total 7*7* cross section [l3 - [171 ]. On the other hand, the successful description of all inclusive and diffractive deep 
inelastic data from HERA by saturation models [l8l432| suggests that these effects might become important in the 
energy regime probed by current colliders. This motivated the generalization of the saturation model to two-photon 
interactions at high energies performed in Ref. (33[, which has obtained a very good description of the data on the 
77 total cross section, on the photon structure function at low x and on the 7*7* cross section. The formalism used 
in Ref. [33[ is based on the dipole picture with the 7*7* total cross sections being described by the interaction 
of two color dipoles, in which the virtual photons fluctuate into (For previous analysis using the dipole picture see, 
e.g., Refs. [13, HH|)- The main assumption made in [33| is that the dipole-dipole cross section can be expressed in 
terms of dipole-proton cross section with the help of the additive quark model. This is a strong assumption which 
deserves a more detailed analysis. This is our first goal in this paper. In particular, we propose a more sophisticated 
connection between the dipole - dipole and dipole - proton scattering amplitudes. Our second goal is to compare the 
predictions obtained using [33| with those obtained with our approach. We also discuss the dependence of the results 
on the dipole scattering amplitude. Here we make use of the state-of-the art parametrization of the dipole scattering 
amplitude 0] ■ We compare the results with the currently available experimental data and provide estimates of the 
total cross sections and photon structure functions which will be measured in the future linear colliders. It is important 
to emphasize that in [33j the cross sections were estimated considering the GBW model [18j], which is inspired on 
saturation physics, and during the last years an intense activity in the area resulted in more sophisticated dipole - 
proton cross sections 19-32], which could be used to estimate the 77 cross sections. In our study we also compare the 
predictions obtained using the solution of the BK equation with those from the phenomenological saturation model 
proposed in [25[ with free parameters updated in 31]. 

Before introducing the required formulas, some comments are in order. Firstly, in our study the heavy quark 
contribution is not included, since the solution of the BK equation used in our calculations 0] has its free parameters 
fixed disregarding the contribution of heavy quarks to the inclusive and longitudinal structure functions. Recently, 
this solution was improved by including the charm and bottom contributions to these observables, which have strong 
effects on the fit parameters 36]. As this solution is not yet public, we postpone for a future study the discussion of 
heavy quark production in 77 interactions. For consistence we also consider the phenomenological saturation model 
proposed in [251 ] without the inclusion of heavy quarks. However, we also consider the updated version obtained in 
[311 ] . where the free parameters were fixed considering the more recent HI and ZEUS data. Secondly, we are assuming 
in our study that fluctuations and correlations produced by the BK equation in the dilute regime (see, for instance, 
Refs. [37]), where this equation reduces to the BFKL equation, can be neglected in the calculation of the dipole - 
dipole scattering cross section. This is a strong assumption. However, results obtained using a toy model in 38], 
indicate that these effects are tamed by saturation in the high-density regime and by the running of coupling in the 
dilute regime. As in our model the scattering amplitude is given in terms of the solution of the BK equation including 
running coupling corrections, we expect the contribution of fluctuations to be small. 



II. THE DIPOLE PICTURE FOR THE TWO-PHOTON CROSS SECTION 



Let us start presenting a brief review of the two-photon interactions in the dipole picture. At high energies, the 
scattering process can be seen as a succession in time of two factorizable subprocesses (See Fig. [T]): i) the photon 
fluctuates into quark-antiquark pairs (the dipoles), ii) these color dipoles interact and produce the final state. The 
corresponding cross section is given by 

<Tij(W 2 ,QlQl) = / dzi / ^n|*?(«i,ri)| a f dz 2 J dV 2 |v^(z 2 ,r 2 )| 2 <Un,r 2 ,n (1) 

a. 6—1 




FIG. 1: The diagram illustrating the 7*7* interaction in the dipole representation. See formula Q. 



where W 2 is the collision center of mass energy squared, the indices i,j label the polarisation states of the virtual 
photons, i.e. T or L, r denotes the transverse separation between q and q in the color dipole, z is the longitudinal 
momentum fraction of the quark in the photon and Y is the available rapidity interval, Y » h^W^/QiC^)- The wave 
functions \^^(z,r)\ 2 and ^{(z, r)| 2 are given by 

|*/(*,r)| a = ^e}{[z 2 + {l-zf]e)Kl{e f r) +m )Kl{e f r)} , (2) 

\¥ L {z,r)\ 2 = 6 -^e 2 f [4Q 2 z*(l-zfK$(e f r)], (3) 

with (e/) 2 = z(l — z)Q 2 + m 2 , e/ and m/ denote the charge and mass of the quark of flavor / and the functions 
Kq and K\ are the McDonald-Bessel functions. Moreover, cr^ b (rx, rs, Y) are the dipole-dipole total cross-sections 
corresponding to their different flavor content specified by the a and b indices. 

The main input for the calculation of the total cross section is the dipole - dipole cross section, which in the 
approximation of two gluon exchange between the dipoles is given by (See, e.g., appendix A from [391 ]) 

a drf (r 1 ,r 2 ,y) = 2™ 2 r 2 [l+ln(-)] (4) 

r 

where r — Min(n,r2) and R = Max(n , r%) , and is energy independent. In the BFKL approach the dipole - dipole 
cross section reads [l| 



(5) 



where x is the BFKL characteristic function, which satisfies the property x(l) — x(l — l)- It implies that (Jjzfkl IS 
symmetric under the exchange r% < — > r 2 of the two dipoles. 
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The behavior predicted by the BFKL equation implies that the cross section violates the unitarity at high energies. 
Consequently, unitarity corrections should be considered in order to tame the BFKL growth of the dipole scattering 
amplitude. In this problem was addressed considering independent multiple scatterings between the onia within 
the color dipole picture, with unitarization obtained in a symmetric frame, like the center-of-mass frame. As demon- 
strated in [4l|, these unitarity corrections can also be estimated considering the Color Glass Condensate formalism, 
which provides a description of the non-linear effects in the hadron wavefunction. It is important to emphasize that 
in general the applications of the CGC formalism to scattering problems require an asymmetric frame, in which 
the projectile has a simple structure and the evolution occurs in the target wavefunction, as it is the case in deep 
inelastic scattering. Therefore the use of the solution of the BK equation in the calculation of the dipole - dipole 
scattering cross section is not a trivial task. Another aspect that deserves a more attention is related to the impact 
parameter dependence of the scattering amplitude. In the eikonal approximation the dipole - dipole cross section can 
be expressed as follows 



a M (r 1 ,r 2 ,Y)=2 J d 2 bN(r u r 2 ,b,Y) 



(0) 



where Af(r\, r 2 , b, Y) is the scattering amplitude for the two dipoles with transverse sizes r\ and r 2 , relative impact 
parameter b and rapidity separation Y. The scattering amplitude M is related to the ^-matrix by S — 1 — A/", with 
the unitarity of the S- matrix implying Af < 1. This constraint is obeyed by the solution of the BK equation. However, 
the dipole - dipole cross section can still rise indefinitely with the energy, even after the black disk limit (Af = 1) has 
been reached at central impact parameters. It occurs due to the non-locality of the evolution, that keeps expanding 
the gluon distribution in the target towards larger impact parameters. This radial expansion is expected to occur 
logarithmically with the energy, in agreement with the Froissart bound,_which should be contrasted with the power 
like rising towards the blackness at fixed impact parameter. Following [42j |. we will assume that the radial expansion 
only affects the subleading energy dependence of <j dd and study the approach towards unitarity limit at a fixed value 
of the target size. Moreover, we will assume that only the range b < R, where R = Max(ri, r 2 ), contributes for the 
dipole - dipole cross section, i.e. we will assume that Af is negligibly small when the dipoles have no overlap with 
each other (b > R). Therefore, we propose that the dipole-dipole cross section can be expressed as follows 

r R 

a dd (r u r 2 ,Y) = 2N(r,Y) / d 2 b — 2irR 2 N(r, Y) , (7) 

Jo 

where N(r,Y) is independent of the impact parameter and satisfies the unitarity bound. The explicit form of a dd 
reads 

a dd (r 1 ,r 2 ,Y) = 2Trr 2 N(r 2 ,Y 2 )9(r 1 -r 2 ) + 2irr 2 N(r 1 ,Y 1 )9(r 2 -r 1 ) , (8) 
where Y^ — ln(l/xi) and 

_ gf + 4m| 

1 ~ Wi + Q 2 ' [ ' 

In contrast, in the phenomenological model proposed in [33| the dipole-dipole cross section was assumed to be given 
by 

^(n.ra.y) = af N(r u r 2 ,Y) (10) 

with ctq' 6 = (2/3)<7o, where ao is a free parameter in the saturation model considered, fixed by fitting the DIS HERA 
data. This relation can be justified by the quark counting rule, as the ratio between the number of constituent quarks 
in a photon and the corresponding number of constituent quarks in the proton. Moreover, it was assumed that 
N(r u r 2 ,Y) = A(r cff , Y = ln(l/x a6 )), where 



2 _ rjrj _ Qf + Q 2 2 + Amj + 4m b 

^ - rlTr 2 dnd Xab W^ + Q 2 + Q 2 



and x ab = ^z^-r a „r b • (ID 



In what follows we will calculate the total 77, 7*7* cross-sections and the real photon structure function F^(x, Q 2 ) 
considering these two models for the dipole-dipole cross section. Before, in the next Section, we discuss in more detail 
the scattering amplitude N(r, Y) used in our calculations. 
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FIG. 2: Pair separation dependence of the rcBK (solid line) and IIM (dashed line) scattering amplitudes at different values of 
x: (a) x = 10~ 3 and (b) x = 10" 6 . 



III. THE FORWARD DIPOLE SCATTERING AMPLITUDE 



The forward scattering amplitude Af(x, r) is a solution of the Balitsky-Kovchegov (BK) equation, which is given in 
leading order by 



dN(r, Y) 
~8Y 



J dri K LO (r, n , r 2 )[N(n ,Y)+N{r 2 ,Y)~ N(r, Y) - N(n , Y)N(r 2 , Y)] , (12) 



where Y = ln(xo/x) (xq is the value of x where the evolution starts), and r% — r — r\. K ho is the evolution kernel, 
given by 

K h °(r, ri ,r 2 ) = ^^, (13) 

where a s is the strong coupling constant. This equation is a generalization of the linear BFKL equation (which 
corresponds of the first three terms) , with the inclusion of the (non-linear) quadratic term, which damps the indefinite 
growth of the amplitude with energy predicted by BFKL evolution. The leading order BK equation presents some 
difficulties when applied to study DIS small-x data. In particular, some studies concerning this equation (43Tj47j have 
shown that the resulting saturation scale grows much faster with increasing energy (A ~ 0.5 ) than that extracted 
from phenomenology (A ~ 0.2 — 0.3). The calculation of the running coupling corrections to the BK evolution kernel 
was explicitly performed in 0, H[, where the authors included a s Nf corrections to the kernel to all orders. In Q the 
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Model 


N(r,Y) 


m (MeV) 


A (MeV) 


Model 1 


rcBK 


198 






IIM-S 


205 




Model 2 


rcBK 




210 




IIM-S 




230 



TABLE I: Parameters used in the calculations. 



improved BK equation was numerically solved replacing the leading order kernel in Eq. (|12p by the modified kernel 
which includes the running coupling corrections and is given by @ 



K 



Bal, 



r, r 1 ,r 2 ) 



N c a s (r 2 ) 
2tt 2 



1 (a s {r\) 
a s (r|) 



1 { a.(rj) 
a s (rf) 



(14) 



Numerical studies of the improved BK equation [5J have confirmed that the running coupling corrections lead to a 
considerable slow-down of the evolution speed, which implies, for example, a slower growth of the saturation scale 
with energy, in contrast with the faster growth predicted by the LO BK equation. Since the improved BK equation 
has been shown to be quite successful when applied to the description of the ep HERA data for the proton structure 
function, we feel confident to use it in other physical situations such as 77 colisions. In what follows we make use of 
the public-use code available in |48j. 

The running couplmg_Balitsky - Kovchegov (rcBK) predictions will be compared with those from the parametriza- 
tion proposed in Ref. [25j | and updated in [3l|, which was constructed so as to reproduce two limits of the LO BK 
equation analytically under control: the solution of the BFKL equation for small dipole sizes, r <§; 1/Q s (x), and 
the Levin- Tuchin law for larger ones, r 3> 1/Q s {x). In the updated version of this parametrization [31|, the free 
parameters were obtained by fitting more recent HI and ZEUS data. In this parametrization the dipole forward 
scattering amplitude is given by 



N(x,r) = 



2(7,+ 



ln(2/rQ a ) 



1 — exp 



-a ln^ (brQ s ) 



, for rQ s (x) < 2 . 
for rQ s (x) > 2 . 



(15) 



1.0 

15) 



where a and b are determined by continuity conditions at rQ s (x) = 2, j s = 0.6194, n = 9.9, A = 0.2545, Q§ 
GeV 2 , xq = 0.2131 and Wo = 0.7. Hereafter, we shall call the model above IIM-S. The first line from Eq, 
describes the linear regime whereas the second one describes saturation effects. 

In Fig. [2] we compare the pair separation dependence of the rcBK and IIM-S forward dipole scattering amplitudes 
at distinct values of x. The main difference between these models is the rapid onset of saturation predicted by the 
IIM-S model. In comparison, the rcBK solution predicts a smooth growth, with a delayed saturation of the forward 
dipole scattering amplitude. Basically, the asymptotic saturation regime is only observed for very small values of x, 
beyond the kinematical range of HERA. 



IV. RESULTS AND DISCUSSION 



Before presentin g ou r results let us discuss the free parameters in our calculations. In the phenomenological 
model proposed in [331, denoted in what follows Model 1, the parameter which is adjusted in order to describe the 
experimental data of the total 77 cross section is the mass of the light quarks. In order to describe the normalization 
of the data the authors are obliged to assume values which are larger than those obtained in the fit of the F2 HERA 
data using the same phenomenological saturation model. Here we follow the same procedure and adjust the light 
quark mass in the wavefunctions in order to fit the real total cross section when using the rcBK and IIM-S models for 
the dipole scattering amplitude. On the other hand, when using the dipole scattering amplitude from Eq. ((HJ, denoted 
Model 2 hereafter, we constrain the light quark mass to the values obtained in the original fits of the F2 HERA data 
0, Hlj]. However, in Model 2, due to the quadratic dependence on the size of the larger dipole [See Eq. (J7|)], the 
contribution of large values of r\ and r-i is quite significant in the total cross section. In order to keep our calculations 
in the perturbative regime we cut the integration on the pair separation at a maximum value of the order of the inverse 
perturbative QCD energy scale. In other words, we stop the r\ and r-i integrations at a maximum dipole size, which is 
chosen to be r max — , with A a free parameter in the model which is expected to be « Aqc d ■ This parameter will be 
fitted in order to describe the total 77 cross section data at high energies. In principle we could keep A = Aqcd fixed, 
understanding it as a clean-cut frontier between perturbative and non-perturbative physics. We could then compute 
the photon-photon cross sections, compare them with data and observe how well perturbative QCD works in this 
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FIG. 3: The total 77 cross section as a function of the energy W for different models of dipole-dipole cross section and dipole 
scattering amplitude. 

domain. Discrepancies between theory and data would be attributed to non-perturbative contributions. We expect 
these contributions to be larger in the case of real photons, which have a larger average transverse radius. However 
the uncertainties in the value of Aqcd would make this separation between the perturbative and non-perturbative 
regimes less precise. In the present work we adjust the value of A. If the value required to fit the data would be 
much different (e.g. much smaller) from Aqcd this would be an indication that it is not possible to describe the bulk 
of data only with perturbative QCD. Surprisingly the obtained values of A, shown in Table I, are close to the most 
accepted values of Aqcd, suggesting that the physics of high energy photon-photon scattering is to a large extent 
perturbative. 

In Fig. |21 we present a comparison between the predictions obtained using the two models for a dd and N(r, Y) and 
the current experimental data for the total 77 cross section at high energies, where for instance rcBK(l) indicates 
that we are using the rcBK solution for the scattering and the Model 1 for the dipole-dipole cross section and so on. 
It is important to emphasize that differently from [331 ] . the Reggeon contributions are not included in our calculations, 
since our focus is in the energy range W > 50 GcV, where these contributions are very small. The parameters used 
in our calculations are presented in Table |U As in [33| , the description of the experimental data [49[ using the Model 
1 is only possible if we assume larger values of the light quark mass in comparison to those used in the description of 
the F2 data, where m u ^d,s = 140 MeV. On the other hand, the value of A necessary to describe the experimental data 
is almost equal to Aqcd, in agreement with our expectations. This result can be interpreted as an indication that 
Model 2 for the dipole - dipole cross section captures the main features of the interaction. Finally, the predictions 
of the two models for a dd are very similar in the kinematical range of the experimental data, independently of the 
N(r,Y) considered. However, at W > 110 GeV, the predictions are very distinct. In particular, at W = 1000 GeV 
they differ by > 25 %, with Model 2 predicting smaller values for the total 77 cross section. 
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FIG. 4: The total 7*7* cross sections as a function of the variable Y = ln(W 2 /Q1Q2) for different values of Q 2 (Q 2 — Q 2 = Q 



We can also compute the two-photon cross section for the case Q 2 ~ Q\ (with large Q\ 2 ) corresponding to the 
interaction of two (highly) virtual photons and also for the case Q\ 3> Q\ corresponding to probing the structure of 
virtual (<3| > 0) or rea l (Q2 = 0) photon at small values of the Bjorken parameter x = Qi/(2qi<l2) (Q 2 = —qf)- For 
instance, the structure function F^{x, Q 2 ) of the real photon (Q 2 = 0, Q\ = Q 2 ) is related in the following way to the 
7*7 total cross-sections: 



F2(x,Q 2 ) 



4ir 2 t 



[<t t , t (W 2 ,Q 2 ,Q 2 2 



0) 



<jl,t(W 2 ,Q 2 ,Q 2 2 



0)]. 



(16) 



In Fig. 2]we present our predictions for the total 7*7* cross section for different photon virtualities. We assume that 
Qi = Q\ = Q 2 an( i analyse the dependence of the cross section on the variable Y = ln(I / F 2 /QiQ2)- We can see that 
the cross sections increase with Y and decrease with Q 2 . Moreover, similarly to the real case, the main difference 
between the predictions is associated to the choice of a dd , with Model 1 predicting larger values for the cross section 
and a steeper growth in rapidity. This difference increases at larger values of the photon virtuality, being a factor w 8 
for Y = 10 and Q 2 = 20 GeV 2 . The experimental point in the right panel is from the L3 Collaboration [50l ]. 

Finally, in Fig. [5] we present our predictions for the x dependence of the photon structure function F^(x,Q 2 ) for 
different values of the photon virtualities. The basic idea is that the quasi-real photon structure may be probed by 
other photon with a large momentum transfer. We present in the lower right panel our predictions for the virtual 
photon structure function. Although there exist only very few data on this observable, its experimental study is 
feasible in future linear colliders. Our results predict that F^x^Q 2 ) increases at small- a;, similarly to predictions for 
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FIG. 5: The photon structure function F^ix, Q 2 ) as a function of x for different choices of the virtualities Q\ and Q\. 

the proton structure function. The current experimental data [5l| are described quite well. As it was seen previously 
in the case of cx 77 and c 7 * 7 », Model 1 predicts a steeper growth with the energy and, consequently, at smaller values 
of x, with the difference between the models increasing at larger values of Q 2 . An important aspect is that the 
predictions obtained using Model 1 are almost independent of the scattering amplitude used in the calculations. On 
the other hand, in Model 2, the predictions depend more strongly on N(r, Y). This makes the study of F^ix, Q 2 ) an 
important source of information about the QCD dynamics at high energies. 

V. SUMMARY 

In this paper we have estimated the main observables to be studied in 77 collisions in the future linear colliders 
using the color dipole picture. In this approach the main input is the dipole - dipole cross section, which is determined 
by the QCD dynamics. We have discussed this quantity in detail and have introduced a new relation with the dipole 
scattering amplitude. In our calculations we use the state-of-the-art of the non-linear QCD dynamics, with the dipole 
scattering amplitude given by the solution of the running coupling Balitsky-Kovchegov equation which are compared 
with the predictions of the phenomenological saturation model proposed in [25[ and updated in [3lj . Moreover, we 
compare our predictions with those obtained using the phenomenological model for the dipole - dipole cross section 
proposed in [331 ]. We demonstrate that these models are able to describe the scarce experimental data at currently 
available energies. However, they differ largely at higher energies, which implies that future experimental data could 
be used to constrain the QCD dynamics. 
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